library(DESeq2)
library(tidyverse)
library(limma)
library(vegan)

load("RData_files/initial_fullddsdesigncountsVsdcountsWGCNA.RData")

remove batch effect of genotype

vsd2 <- removeBatchEffect(assay(Vsd), batch = Vsd$Genotype)

formal analysis of variance in distance matrix (PERMANOVA)

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkoREVTZXEyKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShsaW1tYSkKbGlicmFyeSh2ZWdhbikKCmxvYWQoIlJEYXRhX2ZpbGVzL2luaXRpYWxfZnVsbGRkc2Rlc2lnbmNvdW50c1ZzZGNvdW50c1dHQ05BLlJEYXRhIikKCgpgYGAKCiMgcmVtb3ZlIGJhdGNoIGVmZmVjdCBvZiBnZW5vdHlwZQoKYGBge3J9CnZzZDIgPC0gcmVtb3ZlQmF0Y2hFZmZlY3QoYXNzYXkoVnNkKSwgYmF0Y2ggPSBWc2QkR2Vub3R5cGUpCmBgYAoKCiMgZm9ybWFsIGFuYWx5c2lzIG9mIHZhcmlhbmNlIGluIGRpc3RhbmNlIG1hdHJpeCAoUEVSTUFOT1ZBKQpgYGB7cn0KYWQ9YWRvbmlzMih0KHZzZDIpflRyZWF0bWVudCxkYXRhPWRlc2lnbixtZXRob2Q9Im1hbmhhdHRhbiIscGVybXV0YXRpb25zPTFlNikKYWQKc3VtbWFyeShhZCkKYXMuZGF0YS5mcmFtZShhZCkgJT4lIAogIHdyaXRlX2NzdigiUEVSTUFOT1ZBX2JhdGNoZWZmZWN0cmVtb3ZlZC5jc3YiKQpgYGAKCg==